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Abstract — In this paper we expand our recently introduced 
concept of UW-OFDM (unique word orthogonal frequency di- 
vision multiplexing). In UW-OFDM the cyclic prefixes (CPs) 
are replaced by deterministic sequences, the so-called unique 
words (UWs). The UWs are generated by appropriately loading 
a set of redundant subcarriers. By that a systematic complex 
number Reed Solomon (RS) code construction is introduced 
in a quite natural way, because an RS code may be defined 
as the set of vectors, for which a block of successive zeros 
occurs in the other domain w.r.t. a discrete Fourier transform. 
(For a fixed block different to zero, i.e., a UW, a coset code 
of an RS code is generated.) A remaining problem in the 
original systematic coded UW-OFDM concept is the fact that 
the redundant subcarrier symbols disproportionately contribute 
to the mean OFDM symbol energy. In this paper we introduce 
the concept of non-systematic coded UW-OFDM, where the 
redundancy is no longer allocated to dedicated subcarriers, but 
distributed over all subcarriers. We derive optimum complex 
valued code generator matrices matched to the BLUE (best 
linear unbiased estimator) and to the LMMSE (linear minimum 
mean square error) data estimator, respectively. With the help of 
simulations we highlight the advantageous spectral properties 
and the superior BER (bit error ratio) performance of non- 
systematic coded UW-OFDM compared to systematic coded 
UW-OFDM as well as to CP-OFDM in AWGN (additive white 
Gaussian noise) and in frequency selective environments. 

Index Terms — Cyclic prefix (CP), Estimation, Minimum mean 
square error (MMSE), OFDM, Unique word OFDM (UW- 
OFDM), Reed Solomon coded OFDM. 

I. Introduction 

IN m, (2) we introduced an OFDM (orthogonal frequency 
division multiplexing) signaling scheme, where the usual 
cyclic prefixes (CPs) Q are replaced by deterministic se- 
quences, that we call unique words (UWs). A related but - 
when regarded in detail - also very different scheme is KSP 
(known symbol padded)-OFDM [4], [5|. Fig.QJ -QJ compare 
the CP-, KSP-, and UW-based OFDM transmit data structures. 

Christian Hofbauer was funded by the European Regional Development 
Fund and the Carinthian Economic Promotion Fund (KWF) under grant 
20214/15935/23108. He is currently funded by the Austrian Science Fund 
(FWF): I683-N13. 

The authors are with the Institute of Networked an Embedded Systems, 
Alpen-Adria-Universitat Klagenfurt, Austria (e-mail: mario.huemer@uni- 
klu.ac.at; chris.hofbauer@uni-klu.ac.at), and with the Institute for Informa- 
tion Transmission, University of Erlangen-Nuremberg, Germany (e-mail: 
huber@lnt.de), respectively. 

Copyright (c) 2011 IEEE. Personal use of this material is permitted. 
Permission from IEEE must be obtained for all other uses, in any current or 
future media, including reprinting/republishing this material for advertising or 
promotional purposes, creating new collective works, for resale or redistribu- 
tion to servers or lists, or reuse of any copyrighted component of this work 
in other works. 

Digital Object Identifier 10. 1109/TSP.201 1.2168522 



In CP- as well as in UW-OFDM the linear convolution of 
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(a) Data structure using CPs 
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Fig. 1. Transmit data structure using a CP (a), a KS (b) or a UW (c). 

the transmit signal with the channel impulse response is 
transformed into a cyclic convolution such that the discrete 
Fourier transform (DFT) can diagonalize the channel in the 
frequency domain. However, different to the CP, the UW is 
part of the DFT-interval as indicated in Fig. Q] Furthermore, 
the CP is a random sequence, whereas the UW is deterministic. 
Hence, the UW can optimally be designed for particular 
needs like synchronization and/or channel estimation purposes 
at the receiver side. The broadly known KSP-OFDM uses 
a structure similar to UW-OFDM, since the known symbol 
(KS) sequence is deterministic as well. The most important 
difference between KSP- and UW-OFDM is the fact, that 
the UW is part of the DFT interval, whereas the KS is not. 
The generation of the UW within the DFT-interval introduces 
correlations among the subcarriers, which can advantageously 
be exploited by the receiver to improve the BER (bit error 
ratio) performance. Whilst in both schemes the deterministic 
sequences can be used for synchronization and channel es- 
timation purposes, these correlations are not present in KSP- 
OFDM. We notice that KSP-OFDM coincides with ZP-OFDM 
(zero padded OFDM) (6), if the KS sequence is set to zero. 

Since UW-OFDM time domain symbols contain a block 
of fixed samples, i.e., the UW, the set of all corresponding 
vectors in discrete frequency domain forms a coset to a Reed 
Solomon code (RS code). Usually RS codes of length n are 
defined for a finite field Fq using an element w £ Fq of 
order n, n ■ I = Q — 1, with n,l,Q £ N to define a discrete 
Fourier transform Fq — > Fq in Fq. The set of codewords is 
specified by the fact, that the (inverse) DFT of all codewords 
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contains a block of d m i n — 1 successive zeros, where d m i n 
is the minimum Hamming distance of the RS code. If this 
block of d m i n — 1 successive symbols differs from zero, but 
is also fixed for all codewords, a coset code to an RS code 
is generated in the other domain w.r.t. this Fourier transform 
with the same minimum distance d m i n , c.f. 0. 

All these definitions apply for the field of complex numbers 
and a usual DFT of length n as well. Thus, in UW-OFDM the 
set of frequency domain data vectors defines a coset code to 
an RS code in a quite natural way. In contrast to the usual 
approach using RS codes over a finite field, e.g. F 2 s, for an 
outer code in a concatenated code scheme, in UW-OFDM we 
have an inner RS code over the field of complex numbers if 
any further channel coding scheme is applied, i.e., the OFDM 
guard space is here additionally exploited for redundancy of 
an inner channel coding scheme in a natural way. In (8) we 
showed that algebraic decoding of the introduced complex 
number RS code leads to solving an ill-conditioned system 
of equations which is extremely sensitive to noise. It turns 
out that the application of data estimation approaches like the 
BLUE (best linear unbiased estimator) or the LMMSE (linear 
minimum mean square error) estimator, cf. |9j], is much more 
appropriate than algebraic decoding. 

For SC/FDE (single carrier/frequency domain equalization) 
systems [ 10]-[24|, the benefits of UW based transmission have 
already sufficiently been studied El-lIlQl, ED-ED . The 
introduction of UWs in SC/FDE systems is straightforward, 
since the data symbols as well as the UW symbols are defined 
in time domain. In UW-OFDM the data symbols are defined 
in frequency domain, whereas the UW symbols are defined 
in time domain, which leads to some difficulties. In ll25ll we 
discussed the similarities and differences of the UW approach 
for OFDM and SC/FDE. 

In our concept described in Q), |9j we suggested to 
generate UW-OFDM symbols by appropriately loading so- 
called redundant subcarriers. The minimization of the energy 
contribution of the redundant subcarriers turned out to be a 
challenge. We solved the problem by generating a zero UW 
in a first step, and by adding the desired UW in a separate 
second step. We showed that this approach generates OFDM 
symbols with much less redundant energy [2| than a single 
step or direct UW generation approach as e.g., described in 
[ 18 1. Additionally, we optimized the positions of the redundant 
subcarriers to further reduce their energy contribution. Several 
other attempts of applying UWs in OFDM systems can be 
found in the literature, e.g. in I26l - ll27ll . In all those approaches 
the guard interval and thus the UW is not part of the DFT- 
interval. Therefore, and in contrast to our UW-OFDM concept 
no coding is introduced by these schemes. 

Our systematic complex number RS coded UW-OFDM 
concept presented in HI and shortly reviewed in Sec. [TT] of 
the present paper still suffers from a disproportionately high 
energy contribution of the redundant subcarriers. In ll28l we 
tackled this problem by increasing the number of redundant 
subcarriers while keeping the length of the UW constant. On 
the one hand this approach in fact leads to a reduction of 
the redundant energy contribution and to an improved BER 



performance, but on the other hand the bandwidth efficiency 
decreases compared to the original concept. In ll29ll we in- 
troduced another approach that also focuses on the redundant 
energy contribution. Here we achieved the reduction of the re- 
dundant energy by allowing some systematic noise within the 
guard interval. This method clearly outperforms the original 
UW-OFDM approach, however, a remaining penalty is the fact 
that the UW is disturbed to some extent. In the present paper 
we introduce a different and much more favorable approach to 
overcome the shortcomings of the original UW-OFDM con- 
cept. We no longer primarily focus on the redundant energy re- 
duction. Instead, we suggest to distribute the redundant energy 
over all subcarriers, and we define cost functions that take the 
overall transceiver performance (including the data estimation) 
into account. The corresponding UW-OFDM symbol genera- 
tion procedure introduces a non-systematic complex number 
RS code construction (cf. Sec. ITITb which can be described by 
appropriate code generator matrices. For the data estimation 
we apply two different approaches, namely the BLUE and 
the LMMSE estimator. Sec. [IV] is dedicated to the solutions 
of the arising optimization problems. At first we solve the 
optimization problems numerically, thereafter we analytically 
derive a number of highly interesting general properties of 
optimum code generator matrices and the implications for the 
overall system approach. Moreover, we discuss the properties 
of two particular numerically found code generator matrices. 
In Sec. [V] we show that non-systematic coded UW-OFDM 
can be converted into a UW-SC/FDE system by choosing a 
specific constructed optimum code generator matrix. Finally, 
simulation results are presented in Sec. [VI] We compare the 
novel UW-OFDM approach against our original systematic 
coded UW-OFDM concept and against a classical CP-OFDM 
system, as a reference system we use the IEEE 802.11a 
WLAN (wireless local area network) standard. The spectral 
advantages are discussed, and BER simulation results are 
presented for the AWGN channel as well as for frequency 
selective indoor scenarios. For the latter case we additionally 
investigate the impact of channel estimation errors on the BER 
performance. The results highlight the advantageous properties 
of the proposed scheme. 



Notation: Lower-case bold face variables (a, b,...) indicate 
vectors, and upper-case bold face variables (A, B,...) indicate 
matrices. To distinguish between time and frequency domain 
variables, we use a tilde to express frequency domain vectors 
and matrices (a, A,...), respectively. We further use R to denote 
the set of real numbers, C to denote the set of complex 
numbers, I to denote the identity matrix, (-) T to denote 
transposition, (•)* to denote complex conjugation, (-) H to 
denote conjugate transposition, E[-\ to denote expectation, 
and tr{ } to denote the trace operator. For all signals and 
systems the usual equivalent complex baseband representation 
is applied. 
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II. Review of Systematic Coded UW-OFDM 

A. Unique Word Generation 

We briefly review our original approach of introducing 
unique words in OFDM time domain symbols, for further 
details see Q], (2). Let x n G C iv " xl be a predefined sequence 
which we call unique word. This unique word shall form the 
tail of each OFDM time domain symbol vector. Hence, an 
OFDM time domain symbol vector, as the result of a length- 
7V-IDFT (inverse DFT), consists of two parts and is of the 
form [xj x^] T G C Nxl , whereat only x rf G C^ N ~ N ^ xl is 
random and affected by the data. In the concept suggested in 
HI, El we generate an OFDM symbol x = [xj T ] T with 
a zero UW in a first step, and we determine the final transmit 
symbol x' = x + [0 T xj] by adding the desired UW in 
time domain in a second step. As in conventional OFDM, 
the QAM data symbols (denoted by the vector d G C NdXl ) 
and the zero subcarriers (at the band edges and at DC) are 
specified as part of the frequency domain vector x, but here 
in addition the zero word is specified in time domain as part of 
the vector x = F 
with elements 

The system of equations x = F^x with the introduced 
features can, e.g., be fulfilled by spending a set of redundant 
subcarriers. We let the redundant subcarrier symbols form the 
vector r G <C NrXl with N r = N u , we further introduce a 
permutation matrix P G (C( Nd+Nr * >x ( Nd+Nr \ and form an 
OFDM symbol (containing N — Nd — N r zero subcarriers) 
in frequency domain by 



-l- 

' N 
[FjvljW 



x. Fat denotes the length-iV-DFT matrix 
= e -i¥ kl for kj = 0,1,..., N - 1. 



x = BP 



(1) 



B G (^Nx(N d +N r ) j nser ts the usual zero subcarriers. It 
consists of zero-rows at the positions of the zero subcarriers, 
and of appropriate unit row vectors at the positions of data 
subcarriers. We will detail the reason for the introduction of 
the permutation matrix P and its specific construction shortly 
below. The time - frequency relation of the OFDM symbol 
(before adding the desired UW) can now be written as 
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(3) 



where My are appropriate sized sub-matrices, it follows that 
M 2 id + M 22 r = 0, and hence r = -M 22 1 M 2 id. With the 
matrix 



M^Mai G 



(4) 



the vector of redundant subcarrier symbols can thus be deter- 
mined by the linear mapping 



Td. 



(5) 



Equation (O introduces correlations in the vector x of fre- 
quency domain samples of an OFDM symbol. The construc- 
tion of T, and thus also the energy of the redundant subcarrier 



symbols, depends on the choice of P. The mean symbol 
energy E x > — E[(x.') H x'} can be calculated to 

1 



^(N d a 2 d +a 2 MTT H )) 



(6) 
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cf. J2). E^/N and Ef/N describe the contributions of the data 
and the redundant subcarrier symbols to the total mean symbol 
energy before the addition of the desired UW, respectively, 
and E Xu describes the contribution of the UW. It turns out 
that the energy contribution Ef/N of the redundant subcarrier 
symbols almost explodes without the use of an appropriate 
permutation matrix, or equivalently for P = I. In |[T1 we 
therefore suggested to choose P by minimizing the symbol 
energy E x > or equivalently by minimizing the energy-based 
cost function 



JK = ^tr{TT"} 



(7) 



Note that T is derived from © and ©. 

Example 1: For the parameter choice N = 64, N u = 16, 
and an index set of the zero subcarriers given by {0, 27, 
28, ...,37} (these parameters are taken from the IEEE 802.11a 
WLAN standard ED, see also Table |B, we have N r = 16 
and Nd — 36. The optimum index set for the redundant 
subcarriers as a result of minimizing the cost function in (0 
is {2, 6, 10, 14, 17, 21, 24, 26, 38, 40, 43, 47, 50, 54, 58, 
62}; cf. (TJ. This choice can easily also be described by (Q]) 
with appropriately constructed matrices B and P. We assume 
uncorrelated and zero mean QAM data symbols with the 
covariance matrix C^j = a 2 d \. Fig. |2] shows the mean power 
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Fig. 2. Mean power of individual subcarrier symbols for Example 1, 

values of all individual subcarrier symbols for the chosen 
parameter setup for the case the UW is the zero word x lt = 
and for cr? = 1. The optimized mean power values of the 
redundant subcarrier symbols are the elements of the vector 
cr^diag (TT H ) evaluated for the optimum permutation matrix 



B. Interpretation as a Systematic Complex Valued Reed- 
Solomon Code 



With 



G = P 



£(N d +N r )xN d 



(8) 
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Fig. 3. Codeword generator for the systematic code described by G. 



we can interpret 



Gd 



(9) 



(c e <c( N d+N r )xi^ as a coc [ eworc i G f a systematic complex 
number Reed Solomon code construction with the code gen- 
erator matrix G. As already mentioned above an RS code 
with minimum Hamming distance d m in may be defined as 
the set of codewords, which all show a block of d m in — 1 
consecutive zeros in their spectral transform w.r.t. a Fourier 
transform defined in the (elsewhere usually finite) field from 
which the code symbols are taken; c.f. 0. Here, simply 
time and frequency domains are interchanged and the field 
is the set of complex numbers. Fig. [3] graphically illustrates 
the generation of a codeword c = [cq, ci, CN d +N r -i] T ■ 
Using (O, and with the frequency domain version of the UW 
x lt = F^v [0 t x^] the transmit symbol can now also be 
written as 

x' =F iv 1 (BGd + x u ). (10) 

C. System Model and Preparatory Steps 

After the transmission over a dispersive (e.g., multipath) 
channel a received OFDM time domain symbol can be mod- 
eled as 



H c x' + n 

H c F^ 1 (BGd + x u ) + n, 



(11) 
(12) 



cf. ( TTOb . where n £ C Nxl represents a zero-mean Gaussian 
(time domain) noise vector with the covariance matrix cr^I, 
and H c e C WxA ' denotes a cyclic convolution matrix originat- 
ing from the zero-padded vector of channel impulse response 
coefficients h c G C Nxl . After applying a DFT to obtain 
y r = Fjvy r , we exclude the zero subcarriers from further 
operation, which leads to the down-sized vector y^ = B T y r 
withy d eC(^ +Ar "' xl : 



\ n, r N (BGd - 



x„) + B T F w n. 



(13) 



y<2 = B F/vH r F 

The matrix H c = FatHcF^ 1 is diagonal and contains the 
sampled channel frequency response on its main diagonal. 
H 



B^H^B with H e C( N *+ N r)*(N d +N r ) is a 



down-sized version of the latter excluding the entries corre- 
sponding to the zero subcarriers. The received symbol can now 
be written in the form of the affine model 



Yd = HGd 



HB T x,, 



B T F Ar n. 



(14) 



Note that (assuming that the channel matrix H or at least 
an estimate of the same is available) HB T x„ represents the 
known portion contained in the received vector yd originating 
from the UW. As a preparatory step to the data estimation pro- 
cedure the UW influence is subtracted to obtain the corrected 
symbol y = yd — HB T x u in the form of the linear model 



HGd 



(15) 



with the noise vector v = B T FArn. The vector y serves as 
the input for the data estimation (or equalization) procedure. 
In the following we will consider linear data estimators of the 
form ^ 

5 = Ey, (16) 

where E £ C NdX< ' Nd+N '^ describes the equalizer. 

D. Optimum Linear Data Estimators 

One way to look for an optimum data estimator is to 
assume the data vector to be deterministic but unknown, and 
to search for unbiased estimators. In order for the estimator 
to be unbiased we require 

E[d] = E[Ey] = E£[HGd + v] = EHGd = d. (17) 

Consequently, the unbiased constraint takes on the form 

EHG = I. (18) 

Equ. ( TT8l is equivalent to the ZF (zero forcing) criterion for 
linear equalizers. The optimum solution which is commonly 
known as the best linear unbiased estimator, and which is 
equivalent to the optimum linear ZF equalizer, is found by 
applying the Gauss-Markov theorem, cf. [31], to the linear 
model in HSi . The solution is given by 



Eblue - (G /f H if C^ 1 HG)- 1 G if H^Cr 1 . 



(19) 



We note that since the noise in ( fT5l ) is assumed to be Gaussian, 
d~T9b is also the MVU (minimum variance unbiased) estimator. 
With the noise covariance matrix Cyy = E [vv ff ] = Na^I 
we obtain 



E 



BLUE 



= (G ff H ff HG)-'G il H 



(20) 



The covariance matrix of d = EBLUEy\ or equivalently, the 
covariance matrix of the error e = d — d immediately follows 
to 

Cgg = ^(G^H^HG)- 1 . (21) 

The most common linear data estimator is the LMMSE esti- 
mator which belongs to the class of the Bayesian estimators. 
In the Baysian approach the data vector is assumed to be the 
realization of a random vector instead of being deterministic 
but unknown as assumed above. By applying the Bayesian 
Gauss-Markov theorem [31] to (fT5l l. where we now assume d 
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to be the realization of a random vector, the LMMSE equalizer 
follows to 

Elmmse - (G ff H ff HG + ^!l)-i G ff H ff . (22) 

Expression ( f22b shows huge similarity to the BLUE in d20l ). 
For a\ = the LMMSE equalizer and the BLUE are identical. 

The covariance matrix of the error e = d — d is given by 

Ce^iVd^G^HG + ^I)- 1 . (23) 

III. UW Generation by Optimum Non-Systematic 
Coding 

In Sec. [II] we chose the positions of the redundant subcar- 
riers (represented by the choice of the permutation matrix P) 
such that the redundant energy becomes minimum. For that 
we had to minimize the cost function Je in (0. Nevertheless, 
the mean power of the redundant subcarrier symbols is still 
considerably higher than that of the data symbols; cf. Fig. |2] 
In this section we present a novel and completely different 
approach to optimize the overall system performance by 
adapting our original concept as follows: 

1) We give up the idea of dedicated redundant subcarriers, 
and we allow to spread the redundant energy over all 
codeword symbols. 

2) We define new cost functions that additionally take the 
receiver processing into account. Instead of purely fo- 
cusing on the redundant energy, we define performance 
measures based on the sum of the error variances at the 
output of the data estimator. 

It will turn out that this approach significantly outperforms CP- 
OFDM and also our original systematic coded UW-OFDM. 

A. The Idea of Non-Systematic Coding in UW-OFDM 

With the introduction and optimization of the permutation 
matrix P we minimized the energy contribution of the re- 
dundant subcarrier symbols. From the optimum choice of the 
permutation matrix P we learned that the redundant subcarrier 
symbols shall be distributed approximately equidistantly over 
the codeword c; cf. Fig [2] This means that the redundant 
energy is not concentrated in bundles of subcarrier symbols, 
but it is spread out over the codeword. Nevertheless, the 
portions of the redundant energy are only concentrated on 
the dedicated redundant subcarrier symbol positions. However, 
from this equidistant distribution of the redundant energy 
one could guess, that it might make sense to distribute the 
redundancy over all subcarrier symbols. If we do so we can no 
longer speak of dedicated redundant subcarriers, since every 
subcarrier will then carry an amount of redundant energy 
instead. We incorporate this idea into our UW-OFDM symbol 
generation process by replacing G as defined in ([8]l by a code 
generator matrix G (of the same size as G) which spreads the 
redundancy over all codeword symbols. The code described by 
6 can then be interpreted as a non-systematic code since the 
original data symbols d will not appear in the codeword 

c = Gd (24) 



any longer. G distributes portions of a single data symbol over 
all (or at least several) codeword symbols, and it additionally 
adds redundancy. Consequently, and analogical to G, G can be 
interpreted as a mixture of a linear dispersive preprocessor (or 
channel-independent precoder, cf. 0,) and a channel coder. 
However, G significantly differs from G in the specific way 
how data and redundancy are spread over the codeword. 

In the following we will formulate optimization criteria 
from which G shall be derived. Following the way we 
optimized the permutation matrix P, we could again think of a 
redundant energy minimization. However, since the redundant 
energy will now be smeared over all subcarrier symbols it 
is not clear how to enforce this. Therefore, we no longer 
primarily focus on the redundant energy reduction, but we 
aim for optimization criteria that take the complete transceiver 
processing into account. 

B. Transceiver Cost Function for the BLUE 

Clearly the linear data estimators in d20l > and d22l) can also 
be used for non-systematic coded UW-OFDM, we only have 
to substitute G by G. We first focus on the BLUE given by 
( f20b . A possible approach to optimize the overall transceiver 
performance is to choose the code generator matrix G such 
that the sum over the error variances after the data estimation 
becomes minimum. With d2Tb this would lead to the cost 
function 

J = tr{C gg } = Na 2 n tx {(G^HG)- 1 }. (25) 

We are aiming for a code generator matrix design which shall 
be done only once during system design. Because of that 
reason the dependence of the cost function on the particular 
channel H is inappropriate. We therefore suggest to look for 
an optimum G for the case H = I, that is the AWGN channel 
case. J then reduces to 

J = iVCT J 2 l tr{(G ff G)- 1 }. (26) 

In the simulation section we will demonstrate that the finally 
derived non-systematic coded UW-OFDM systems not only 
perform superior in the AWGN channel case, but also and 
particularly in frequency selective channelo As we will 
see this comes from the advantageous combination of the 
channel coding and dispersive preprocessing abilities of the 
optimized code generator matrices. But let us come back to 
the formulation of an appropriate optimization criterion and to 
J as given in (f26j): We could now try to minimize J for a given 
cr 2 , where the particular choice of a 2 is obviously irrelevant. 
However, different choices of 6 lead to different mean OFDM 
symbol energies and consequently to different ratios E s /a^, 
where E s denotes the mean energy of an individual QAM data 
symbol. Since it is not desirable to reach the goal of a small 
sum of error variances at the cost of a huge transmit energy, it 
is much more reasonable and fair to fix the ratio c = E s /a^ 
during the optimization. To obtain an expression for E s we 

'it will turn out, that the solution to the formulated optimization problem 
is not unambiguous. 
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calculate the mean OFDM symbol energy E x (for the case of 
a zero UW) first: 



E x = E[x H x] 

= £[(F^BGd) H F^BGd] 

= ^[d K G ir B T F JV F^ 1 BGd] 

= — E\d H G H Gd] 

N 1 1 

With a H b = tr{ba ff } we can further write 



(27) 



E x = -E[tv{G H Gdd H }) 

= ltr{£[G"Gdd*]} 

- ^tr{G H G}. (28) 

The mean QAM data symbol energy E s follows to E s = 
Ex./N d . With c = E s /al we obtain 



C. Transceiver Cost Function for the LMMSE Data Estimator 

In this subsection we assume that the LMMSE estimator 
will be used for data estimation. Again we are aiming for 
optimizing the overall system performance by minimizing the 
sum over the error variances after the data estimation. Using 
(l23b with H = I this leads to 

J = iV(72tr|(G ff G + ^I)- 1 |. (34) 

Like in the considerations for the BLUE we fix c = E s /o-^, 
and we therefore insert ( f29b into 04b to obtain the cost 
function 



Ji 



LMMSE 



Alternatively, the cost function can also be written as 

\ -l " 



Jlmmse — cr d tr 



G H G + 1 



E s _ ^tr{G ff G} 



tr{G H G} 



(35) 



(36) 



Inserting (1291 into 



c cNN d 

finally yields the cost function 



(29) 



jBLUE = | tr {G«G}tr{(G ff G)-} (30) G = argmin{ J LMMSE } s.t. F^BG 



The constrained optimization problem to find G can finally 
be written as 

o • 



(37) 



The cost function Jblue measures the overall transceiver 
performance at a fixed ratio Eg/cr^. However, the particular 
choice of the ratio c = E s /a^ is obviously irrelevant for the 
searching of optimum code generator matrices. 

An optimum code generator matrix G shall consequently 
be found by minimizing Jblue, but in addition, in order that 
the zero UW is generated in the time domain, G has to be 
constrained to fulfill 



F^BGd 



for every possible data vector d, or equivalently 



F- X BG 



(31) 



(32) 



Hence, every column vector of G has to be orthogonal to 
the N u lowermost row vectors of F^B G C Nx ^ Nd+Nr K 
Note that F^ V 1 B is composed of those columns of F^ 1 that 
correspond to the non-zero entries of the OFDM frequency 
domain symbol x. 

Consequently, we have to solve a constrained optimization 
problem for 6, which can finally be written as 



G = argmin {Jblue } s.t. F^BG 



For sufficiently large c we have Jlmmse ~ Jblue, and the 
particular choice of c is again irrelevant for the searching 
of optimum code generator matrices. However, this is not 
immediately apparent for small values of c. The solutions to 
the optimization problem will lead to code generator matrices 
G matched to the LMMSE 'decoding' procedure. 

IV. Solutions of the Optimization Problems 

In this section we at first solve the optimization problems 
in d33l and d37l > numerically. The solutions are ambiguous, 
however all found code generator matrices share a number of 
common properties which we will discuss in detail. With the 
help of analytical considerations which are partly shifted to the 
appendices we show that all found code generator matrices not 
only correspond to local minima but to the global minimum 
of the associated cost function. 



A. Preparatory Steps and Numerical Solution with the Steep- 
est Descent Algorithm 

In this section we use the steepest descent algorithm to nu- 
merically solve the optimization problems in d3~3l ) and ( f3Tb . As 
a preparatory step we transform the constrained optimization 
problems into unconstrained problems. For that and according 
to ([H) we write G in the form 



(33) 



G = AP 



The solutions to the optimization problem will lead to code 
generator matrices G matched to the BLUE 'decoding' pro- 
cedure. 



(38) 



with a non-singular real matrix A e ^{N d +N r )x{N d +N r ) ^ 
and with a fixed permutation matrix P as e.g., found by 
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minimizing the cost function 


in 


0. 


The constraint in <j33j 


and d37l) can thus be rewritten 


as 








F^BAP 


I 

f 




* 



(39) 



With the introduction of 



2) Random Initialization: In the second approach we chose 
each element of A^ ' as a realization of a Gaussian random 
variable with mean zero and variance one: 



[A<°>]y ~JV(0,1) 

We denote the resulting code generator matrix with G'' 



(46) 



M 

the constraint in 
f as 



F^BAP 



Mil 
M 21 



Mi 2 
M 22 



(40) 



can now simply be fulfilled by choosing 



f =-(M 22 )- 1 M 



21- 



(41) 



For both cost functions Jblue and Jlmmse we can claim the 
following: By using the initialization as in (l44t the steepest 
descent algorithm converges at least one order of magnitude 
faster compared to the case when d46l ) is used. For the random 
initialization approach the resulting code generator matrix 
generally varies from trial to trial. 



That means that for a given non-singular real matrix A, the 
matrices T and G can un-ambiguously be calculated by (|4TT > 
and ( 1381 ), respectively, such that the constraint in d33l and 
d37| i is automatically fulfilled. We can therefore consider the 
cost functions Jblue an d Jlmmse as functions of the real 
matrix A, where T and G have to be determined by (PiTT i and 
(l38l . respectively. The steepest descent algorithm can then be 
applied to the unconstrained optimization problems 



L opt 



argmin {Jblue, lmmse} 



(42) 



in a straight forward manner. In d42l either Jblue or Jlmmse 
is minimized. By following this approach the steepest descent 
algorithm automatically only searches within a subset of 
matrices G that fulfill the constraint in ( |33l and d37| i, or 
in other words that produce a zero UW in the OFDM time 
domain symbols. 

For the steepest descent algorithm the gradients of the 
cost functions Jblue and Jlmmse with respect to the real 
matrix A are required. We approximated the partial derivations 

dJ/d[A] l3 by 



dJ ^ J([A] l3 + e) - J([A] l3 - e) 



(43) 



with a very small e. For J we inserted Jblue or Jlmmse, 
respectively. We used two different approaches for the initial- 
ization of the steepest descent algorithm: 

1 ) Initialization with the Code Generator Matrix G: In our 
first approach we chose the initialization 



A<°> = I 
which implies f (0) = T and 

G (°)=P[I T T ] T = G. 



(44) 



(45) 



The iterative optimization process consequently starts with the 
code generator matrix G of our original systematic coded UW- 
OFDM concept, which can definitely be assumed to be a good 
initial guess. We denote the resulting optimum code generator 
matrix (found after convergence of the algorithm) with G'. 



B. General Properties of Optimum Code Generator Matrices 
for the BLUE 

Interestingly, all found local minima for the BLUE based 
numerical optimization feature the same value of the cost func- 
tion Jblue.iiuih independently of the choice of the initialization 
A' ). Another highly interesting finding is, that all resulting 
code generator matrices (again independently of A^ ') feature 
the property 

G H G = al (47) 

with some constant a (which may vary dependent on the 
results of the optimization process). This property has a 
number of important implications. First, inserting ( f47T > into 
the cost function ( 130b leads to 



J B LUE,min = — J" (N d Ot) (iVrftt" 1 ) 

cN d 



(48) 



which is in agreement with the numerically found local 
minima. We can conclude that every 6 fulfilling 



G ff G = al 



and F^BG = 



(49) 



for any value of a will also result in the same value JBLUE.min 
of the cost function, and will produce a zero UW in time 
domain. Second, if we apply a code generator matrix satisfying 
d49l ), then the error covariance matrix after the data estimation 
in the AWGN channel is given by 



G„„ 



ee.BLUE — 



(50) 



This simply follows from inserting d47| i and d29l ) into ( TSTI ). 
As an important consequence we can conclude that the noise 
at the output of the BLUE is uncorrelated under AWGN 
conditions. This is clearly in contrast to systematic coded 
UW-OFDM, where Cgg is non-diagonal also in the AWGN 
channel case. And third, d47b implies that all singular values 
of 6 are identical. To show this we consider a singular value 
decomposition (SVD) of 6 as 



G = USV 



H 



(51) 



with unitary matrices U and V, and with the matrix £ = 
[D 0] , where D is a real diagonal matrix having the 



x 



singular values si, S2, sjsr d of G at its main diagonal. With 
(|47l i we therefore have 

al = G H G = VE H XJ H VEV H = VD 2 V H 
^aI = D 2 = diag {«?,«!,... ,«^}. (52) 

From d52l it follows that G ff G = al implies a = s 2 = 



*2 — ' " ' ^ 
be written as 



:= s 2 . The property in (|47| > can therefore also 



G H G = s 2 I. 



(53) 



The argumentation can also be done the other way round: If all 
singular values of G are identical then we have G H G = al 
with a = s 2 . 

An open question is still whether the value of the cost 
function JBLUE.min as in (l48l corresponding to the numer- 
ically found local minima depicts the global minimum of 
the constrained optimization problem in (l33l l. To answer this 
question we now at first merely concentrate on the cost 
function Jblue> and we disregard the constraint in d33l for 
a moment: Let s = [si S2 • • • sjv d ] be the vector of 
singular values of G. In Appendix [A] we will analytically 
show that 9Jblue/<9s = if and only if all singular values 
of G are identical. Consequently, every possible candidate G 
for a local minimum satisfies G^G = s 2 I (cf. (1521 and its 
implications). Inserting G H G = s 2 I into the cost function 



leads to the same expression as in 



and hence, every 



G fulfilling G H G = s 2 I results in the same (and minimum) 
value JBLUE,min = & d N d /c which therefore constitutes the 
global minimum of the cost function. 

We now come back to the constrained problem in d33l i: 
From our numerical solutions we know that matrices exist, 
that firstly satisfy G = s 2 I and therefore result in 
the global minimum of the cost function Jblue, and that 



secondly fulfill the constraint F N BG 



With these 



considerations we finally end up with the following important 
proposition: 

Properties of optimum code generator matrices: A code 
generator matrix G is optimum, i.e., leads to a global mini- 



mum of the constrained optimization problem in 
only if G satisfies 

G H G = 



s 2 I 



and 



F^BG = 



;/ and 

(54) 
(55) 



where s := Si = S2 = • • • = Spf d are the (all identical) 
singular values of G. The global minimum of the cost function 
is given by d481 >, and the error covariance matrix after data 
estimation (in the AWGN channel) is the scaled identity matrix 
as given in d501 >. 

Note that because of d54T > the colums of any optimum 
code generator matrix G form an orthogonal basis of an N d - 
dimensional subspace of C^ Ard+A ' r ^ xl . Furthermore, as already 
discussed above, (EBT l implies that every column vector of an 
optimum G is orthogonal to the N u lowermost row vectors 
of F- X B. 



C. General Properties of Optimium Code Generator Matrices 
for the LMMSE Estimator 

In Appendix [B] we will analytically show that 
c?Jlmmse/<9s = if and only if all singular values of G are 
identical. All other findings from Sec. IIV-BI also hold for the 
LMMSE estimator based transceiver optimization, except the 
particular expressions for J L MMSE,min and C^lmmse differ 
slightly. With (07]i, (|23]l and (|36]l it immediately follows that 



J] 



LMMSE, min 



y ee, LMMSE 



c+l' 



1 



(56) 



(57) 



As a consequence of the above findings we learn that a 
code generator matrix which is optimum for the BLUE based 
'decoding' procedure is automatically also optimum for the 
LMMSE based data estimation (and vice versa). 

D. Normalized Optimum Code Generator Matrices 

From (l48T > and ( l56i l we learn that the particular value of 
a = s 2 does not play any important role. In our simulations, 
cf. Sec. IVII we therefore normalized all found code generator 
matrices such that a — s 2 — 1 or G H G = I. The 
columns of any normalized optimum code generator matrix G 
form an orthonormal basis of an A^^-dimensional subspace of 
£^{N d +N r )xi _ ^ s ano th er consequence of s = 1 the operation 
c = Gd becomes energy-invariant and we have Ez = = 
N^a-j, and the mean energy of an OFDM time domain symbol 
(for the zero UW case) follows to E x = N d a 2 d /N; cf. 



E. Comparison of Generator Matrices obtained from different 
Initialization Strategies 

We will now discuss some further interesting properties 
of two particular numerically found solutions. Here we only 
concentrate on G' and on one particular solution for G" (and 
the corresponding matrices A' opt and A" t ) found by applying 
the LMMSE estimator based optimization with c = 1 for the 
initializations as described in (l44l and (|4oT >. respectively. 

The matrix A' opt features the symmetry property 

A o P t = I a o • ' - a k/2-i flip{ a k/2-i} ' • - m P{ a o}]> (58) 



where N a = N d + N r and the a- with i = 0, 1, N a /2 - 1 
are the first N a /2 colums of A' t . The corresponding code 
generator matrix G' shows the same symmetry property as G 
in Section namely 

&=]&■■■ SN d /2-l flip{(g^/ 2 -l)*} ■ • - fli P{(g0)*}]- 

(59) 

Here the g[ with i = 0, 1, Nd/2 - 1 are the first N d /2 
columns of G'. These symmetry properties do not hold 
for A" pt and G", respectively. It appears that the matrices 
A' opt and A" pt show a completely different construction. A^ pt 
approximately features a band matrix structure, where all 
dominant entries are positioned on the main diagonal and 
the first few diagonals directly above and below the main 
diagonal. The remaining elements are close to zero. In contrast 
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A" t is a full matrix. This results in completely different Gsc fulfills 
(pre)coding properties of G' and G". G' can be regarded 
as the natural perfecting of G and is constructed such that the 
energy of one data symbol is mainly (however, not exclusively) 
spread locally. In contrast G" spreads the energy of each data 
symbol approximately uniformly over the codeword c. While 
G' and G" perform identically in an AWGN environment the 
discussed differences lead to a quite different behavior of G' 
and G" in frequency selective environments. We will exem- 
plify these differences and the consequences in the simulation 
section. 

Example 2: We apply the system parameters as in Example 1 
and investigate the mean power levels for all individual subcar- 
riers for the case the non-systematic code generator matrix G' 
is applied. The mean power values for the codeword symbols 
correspond to the diagonal elements of the covariance matrix 
Cgg = Fig. |4] shows the power distribution over 

all subcarrier symbols (additionally also including the zero 
subcarrier symbols) again for the case the UW is the zero 
word x u = and for a\ = 1. We can clearly identify 



Gf^Gsc 



ATI 



and F^BGsc 



(61) 



and consequently constitutes an optimum code generator ma- 
trix; cf. (l54l and ( f55l >. Furthermore, it is apparent that Gsc 
generates a UW-SC/FDE signal with a zero UW x u = since 
the time domain symbol vector follows to 



x = F^BGd 



T 




d 





d = 










(62) 
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Fig. 4. Mean power of individual subcarrier symbols for Example 2. 

that our chosen optimality criterion also implies a significant 
reduction of the power levels of the former redundant sub- 
carriers compared to the original UW-OFDM approach; cf. 
Fig. |2 Furthermore, it can be seen that the redundant energy 
is now smeared over all subcarriers; consequently we can 
no longer speak of data or redundant subcarriers. The power 
levels are quite similar for all symbols, the only exceptions 
are the two subcarrier symbols at the band edges. Due to 
the normalization of G' such that (G') H G' = I we have 
tr{C ££ } = ajtr{G'(G') H } = o\Xx{(G') H G'} = N d a% In 
Fig. [4] the sum over all mean power levels is therefore 36. 



V. On the Relationship between UW-OFDM and 
UW-SC/FDE 

With the help of a specific constructed code generator matrix 
non-systematic coded UW-OFDM can be converted into a 
UW-SC/FDE system. For that we assume for the moment 
that no zero subcarriers are used, or B = I 6 M. NxN and 
Nd = N — N r , and we consider the matrix 



Gsc — F 



A' 



(60) 



Simulation results for a UW-SCFDE system in comparison 
with systematic coded UW-OFDM can be found in |25|. 

VI. Simulation Results 

In this section we present a number of simulation results 
to show the advantageous features of the developed non- 
systematic coded UW-OFDM concept. In our simulations the 
transmitter processing starts with optional (outer) channel 
coding, interleaving and QAM-mapping (we apply QPSK 
symbols unless specified otherwise). We used the same outer 
convolutional encoder as defined in ll30l . and we show results 
for (outer) coding rates r — 3/4 and r = 1/2, respectively. 
Next the complex codewords are determined by either using 
G, G' or G". Here, G' and G" are the particular code 
generator matrices dicussed in Sec. IIV-EI We note that the 
(optional) outer convolutional code is a binary code while the 
inner code described by G, G' or G" is an RS code over 
the field of complex numbers. The latter is naturally always 
inherently present due to the proposed way of generating UW- 
OFDM symbols with zero UWs at their tails. After applying 
a code generator matrix, zero subcarriers are filled in, and the 
IFFT (inverse fast Fourier transform) is performed. Finally, the 
desired UW is added in time domain. At the receiver side the 
processing for one OFDM symbol starts with an FFT, then 
the influence of the UW (HB T xJ is subtracted; cf. (fl4] i- 
(fl5l l. Next the data estimation is applied. Finally demapping, 
deinterleaving and (outer) channel decoding are performed. 
For the applied soft decision Viterbi channel decoder the 
main diagonal of the appropriate matrix Cgg is used to 
specify the (in case of transmitting over frequency selective 
channels) varying noise variances along the subcarriers after 
data estimation. 



A. Simulation Setup 

We compare our UW-OFDM approaches with the classical 
CP-OFDM concept. The IEEE 802.11a WLAN standard |30) 
serves as reference system. We apply the same parameters for 
UW-OFDM as in ll30l wherever possible, the most important 
parameters used in our simulations are specified in Table [I] 
The sampling frequency has been chosen to be f s = 20MHz. 
As in [30 1 the indices of the zero subcarriers within an OFDM 
symbol x are set to {0, 27, 28, ...,37}. The indices of the 
redundant subcarriers are chosen to be {2, 6, 10, 14, 17, 21, 
24, 26, 38, 40, 43, 47, 50, 54, 58, 62} as already discussed 
in the Example 1 in Sec. [ID Note that in conventional CP- 
OFDM like in the WLAN standard, the total length of an 
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TABLE I 

Main PHY parameters of the investigated systems. 





802.11a 


UW-OFDM 


Modulation schemes 


QPSK, 16QAM 


QPSK, 16QAM 


Coding rates (outer code) 


1/2, 3/4 


1/2, 3/4 


Occupied subcarriers 


52 


52 


Data subcarriers 


48 


36 


Additional subcarriers 


4 (pilots) 


16 (redundant) 


DFT period 


3.2 fj,s 


3.2 fis 


Guard duration 


800 ns (CP) 


800 ns (UW) 


Total OFDM symbol duration 


4 /LtS 


3.2 fis 


Subcarrier spacing 


312.5 kHz 


312.5 kHz 


Welch Power Spectral Density Estimate 




Frequency (MHz) 

Fig. 5. Power spectral density comparison: CP-OFDM, UW-OFDM using 
G (systematic coded), G' (non-systematic coded - case 1) and 6" (non- 
systematic coded - case 2). 



OFDM symbol is given by Tuft + Tqi- However, the guard 
interval is part of the DFT period in the UW-OFDM approach 
which leads to significantly shorter total symbol durations. It 
is therefore important to mention that the compared systems 
show (almost) identical bandwidth efficiencies. 

Note that in the IEEE 802.11a standard 4 pilot subcarriers 
are specified. Those are used for estimation and synchro- 
nization purposes at the receiver side. In our UW-OFDM 
approaches we omitted these pilots, because the unique word, 
which is deterministic, shall (at least) take over the estimation 
and synchronization tasks which are normally performed with 
the help of the 4 pilot subcarriers. In order to make a fair 
BER performance comparison, the energy of the UW related 
to the total mean energy of a transmit symbol is set to 4/52 
in our BER simulations. This exactly corresponds to the total 
energy of the 4 pilots related to the total mean energy of a 
transmit symbol in the IEEE standard. As UW we applied a 
linear chirp sequence exhibiting the same bandwidth as the 
data signal, and featuring a constant envelope in time domain 
and approximately a constant envelope in frequency domain. 
However, the particular shape of the UW has no impact on 
the BER behavior; cf. 0. 



B. Power Spectral Density 

Fig- H shows the estimated power spectral densities (PSDs) 
of simulated UW-OFDM bursts and of a CP-OFDM burst. 
For all cases we simulated a burst composed of a preamble 
(in all cases the IEEE 802.11a preamble), and a data part 
consisting of 1000 bytes of data. We used an outer channel 
code with coding rate r = 1/2. For the UW-OFDM concepts 
(G, G', G") we exceptionally applied the zero UW for 
these PSD investigations. Note that we did not use any 
additional filters for spectral shaping. Fig.|5]clearly shows that 
the UW-OFDM spectra feature a significantly better sidelobe 
suppression compared to the CP -OFDM spectrum. The out-of- 
band emissions generated by G and G' are more than 15dB 
below the emissions of the CP-OFDM system. The emissions 
are even notably lower for G". Furthermore, the spectra for 
G' and G" feature an extremely flat in-band region compared 
to systematic coded UW-OFDM. This can be explained by 
the fact, that for the systematic coded case the mean power 
strongly varies between data and redundant subcarriers, cf. Fig. 
while all subcarriers (except the ones at the band edges) 
show almost equal power in the non-systematic case. 

C. BER Simulation Results with Perfect Channel Knowledge 

We will now show BER simulation results for the AWGN 
channel as well as for frequency selective environments. To 
avoid confusions in the figures we at first only use 6' for the 
non-systematic coded UW-OFDM system, and at the end of 
the section we will then also show and interpret results for 
G". Perfect channel knowledge at the receiver is assumed in 
all simulations. 

1) Results for the AWGN case: Clearly, OFDM is designed 
for data transmission in frequency selective environments. 
Nevertheless, we start our comparison with simulation re- 
sults in the AWGN channel, since we optimized the non- 
systematic code generator matrices for that case. In Fig. [6] 
the BER behavior of the IEEE 802.11a CP-OFDM based 
standard, and of both, the systematic coded (G) and the non- 
systematic coded (6') UW-OFDM approach are compared 
under AWGN conditions. No outer code is used for these 
simulations. Simulation results are provided for the BLUE 
and for the LMMSE data estimator, respectively. For the sys- 
tematic coded system additional results for a simple channel 



inversion (CI) receiver (E = 



I 



H x ) are included for 



comparison reasons. For the systematic coded UW-OFDM 
system it can be observed that the BLUE and the LMMSE 
estimator clearly outperform the CI receiver, and the LMMSE 
estimator performs slightly better than the BLUE. For non- 
systematic coded UW-OFDM we only plotted one single curve 
since the BLUE and the LMMSE data estimator perform 
completely identical. This is somewhat surprising since the 
error variances after data estimation are not identical; cf. ( l50l ) 
and (l57l >. However, we found that in the AWGN case the 
QPSK symbol estimates of the BLUE and of the LMMSE 
data estimator always lie in the same decision region of the 
constellation diagram, and the difference in the error variances 
does not translate into a difference in the BER performance. To 
give some numbers we compare the performances at a bit error 
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Fig. 6. BER comparison between the UW-OFDM approaches (G and &') 
and the CP-OFDM system in the AWGN channel (QPSK). 



Fig. 7. BER comparison between the UW-OFDM approaches (G and 6') 
and the CP-OFDM system in a frequency selective environment w/o outer 
coding (QPSK). 



ratio of 10~ 6 . Systematic coded UW-OFDM performs slightly 
worse compared to the CP-OFDM reference system, the non- 
systematic coded UW-OFDM system outperforms CP-OFDM 
by ldB and systematic coded UW-OFDM (with LMMSE 
data estimation) by 1.6dB, respectively. We consider this as 
a remarkable performance of the novel non-systematic coded 
UW-OFDM system. 

2) Results for frequency selective environments (G, G'): 
We now turn to results in frequency selective indoor envi- 
ronments. Since the LMMSE estimator always outperforms 
the BLUE in dispersive channels we only concentrate on the 
LMMSE estimator in the following. For the simulation of 
indoor multipath channels we applied the model described 
in lf32ll . which has also been used during the IEEE 802.11a 
standardization process. The channel impulse responses are 
modeled as tapped delay lines, each tap with uniformly dis- 
tributed phase and Rayleigh distributed magnitude, and with 
power decaying exponentially. The model allows the choice 
of the channel delay spread. For the following simulations we 
have generated and stored 5000 realizations of channel impulse 
responses, all featuring a delay spread of 100ns and a total 
length not exceeding the guard interval duration. Furthermore, 
the channel impulse responses have been normalized such that 
the receive power is independent of the actual channel. The 
subsequent figures represent BER results averaged over that 
5000 channel realizations. 

We start with simulation results for the case no outer code 
is used; cf. Fig. [7] The gain achieved by systematic coded 
UW-OFDM over CP-OFDM is already remarkable. Besides 
the coding gain achieved by G together with the LMMSE 
data estimator this mainly comes from the fact, that due 
to the dispersive preprocessing property of G data symbols 
corresponding to deep fading holes in the channel's frequency 
response can still be detected reasonably, since portions of 
these data symbols are also available at redundant subcarriers. 
Further, the non-systematic coded UW-OFDM system outper- 
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Fig. 8. BER comparison between the UW-OFDM approaches (G and 
6') and the CP-OFDM system in a frequency selective environment with 
additional outer coding (QPSK). 



forms the systematic coded one by another 1.6dB, even though 
6' has been optimized for the AWGN channel case. 

Next we present simulation results for the case the ad- 
ditional outer channel code is used; cf. Fig. [8] For both 
outer coding rates the UW-OFDM approaches outperform CP- 
OFDM, and non-systematic coded UW-OFDM shows by far 
the best performance. The gains of non-systematic coded UW- 
OFDM over CP-OFDM at a bit error ratio of 10~ 6 are 1.9dB 
and 1.7dB for r = § and r = \, respectively, the gains over 
systematic coded UW-OFDM are l.ldB for both coding rates. 

Similar tendencies can also be observed in case 16QAM 
symbols are applied as modulation alphabet; cf. Fig. [9] Non- 
systematic coded UW-OFDM again significantly outperforms 
CP-OFDM by 1.6dB and 1.3dB for r = \ and r = §, 
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Fig. 9. BER comparison between the UW-OFDM approaches (G and 
G') and the CP-OFDM system in a frequency selective environment with 
additional outer coding and 16QAM as modulation alphabet. 

respectively (again measured at a bit error ratio of 10 -6 ). 
However, the gain of systematic coded UW-OFDM over CP- 
OFDM shrinks to 0.2dB for r = h, and even turns to a loss 
of 0.5dB for r = 4. Consequently and quite remarkably, the 
achieved gain of G' over G turns out to be notably larger in 
case of 16QAM compared to QPSK modulation. 

3) Results for frequency selective environments fG', G"): 
Finally we compare the performance of the two different 
derived code generator matrices G' and G". We can state 
that in the AWGN channel they feature exactly the same 
performance. This was expected since every optimum code 
generator matrix shows the same error covariance matrix 
for the AWGN case. However, they feature quite a different 
behavior in dispersive channels. In Sec. IIV-EI we already 
discussed the different structures of G' and the particularly 
chosen G". We remind the reader that G' is constructed 
such that the energy of one data symbol is mainly (however, 
not exclusively) spread locally. G' can be regarded as the 
natural perfecting of G. In contrast 6" spreads the energy of 
each data symbol approximately uniformly over the codeword 
c. From this point of view the system with G" behaves 
comparable to a single carrier system, where the energy of 
each individual QAM symbol is also approximately uniformly 
distributed over the whole bandwidth. In contrast the system 
with G' still rather shows more similarity to classical OFDM, 
where a subcarrier exactly corresponds to one QAM symbol. 
Fig. [10] shows that 6" features extremely good performance 
without an outer code, it significantly outperforms G' in that 
case. For a coding rate of r = |, G" still performs slightly 
better than G', while for r — |, G' clearly outperforms 
G". The coding gain achieved by a strong outer code in a 
frequency selective channel is high for 6' as it might be 
expected for a system rather related to classical OFDM, while 
it is comparably low for G" with its particular dispersive 
preprocessing properties making the transceiver rather behave 




E b /N (dB) 

Fig. 10. Comparison of the different non-systematic coded UW-OFDM 
systems: G' versus G" (QPSK). 

like a single carrier system. 

D. BER Simulation Results for Imperfect Channel Knowledge 

Up till now we have presented performance results assuming 
perfect channel knowledge at the receiver. In this subsec- 
tion we investigate the effect of channel estimation errors 
on the BER performance. Since in UW-OFDM the channel 
H is incorporated into the receiver processing quite in a 
different way as in CP-OFDM, it is not immediately obvious 
whether channel estimation errors will degrade the systems' 
BER performance in the same scale. While in CP-OFDM 
the data estimator is simply given by E = H \ possible 
data estimators for UW-OFDM are given by d20T > and d22b , 
respectively. H is usually replaced by an estimated version 
H, which has a degrading impact on the BER performance. 
In UW-OFDM channel estimation errors have an additional 
impact, namely in the processing step where the influence of 
the UW is subtracted from the received symbol; cf. (IT4l)-(IT5l). 

To investigate the influence of channel estimation errors we 
apply a standard preamble based channel estimation procedure, 
which we briefly describe below. We use the IEEE 802.11a 
preamble defined in [30 1. This preamble contains two identical 
BPSK (binary phase shift keying) modulated OFDM symbols 
(preceded by a guard interval) dedicated to channel estimation 
which we denote by x p = x Pl = x P2 G <C Nxl . Note 
that for the downsized frequency domain version x p .d = 
B t Fatx p we have x p , d e {-1, i}(JVd+JV r )xi_ Let y pi and 
y P2 be the received noisy preamble symbols, and let y p d = 
iB T FAr(y pi +y P2 ). Then a first course unbiased estimate of 
the vector of channel frequency response coefficients naturally 
follows as 

^ fc ]=Ni = y fc ]^w (63) 

for k = 0, (Nd + N r — 1). The latter step is true since 
Xp,d[fc] € { — 1,1}. This course channel estimate can be 
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significantly improved or rather noise reduced by making the 
usually valid assumption that the channel impulse response 
does not exceed the guard duration N u , With the vector of 
channel impulse response coefficients h e C^" x 1 and its zero 



padded version h c £ 



, respectively, this assumption can 



be incorporated by modelling the course channel estimate as 



hi = B T FArh c + n 
h 



B T F 



N 







(64) 



where n G t£(N d +N r )xi re p resen ts a white Gaussian (fre- 
quency domain) noise vector. By decomposing the DFT matrix 



as F N = 



Mi M 2 



with Mi e C 



NxN u 



C Nx(N-N u )^ ^ can be r e W1 -i tt en as 

hi =B T Mih + h. 



and M 2 G 



(65) 



From the linear model in ( I65l l the MVU estimator of the 
channel impulse response follows to 



h = (Mf BB T Mi) 1 Mf Bhi, 



(66) 



cf. IT3T1 . Going back to frequency domain, and again excluding 
the zero subcarriers from further operation, delivers the final 
and highly noise reduced frequency domain channel estimate 



B T F 



N 



= B Mi (Mf BB Mi) MfBhi. (67) 
w 

Note that the smoothing matrix W 6 Nd+N ^ x ( Nd+N ^ 
does not depend on the channel, and has to be calculated only 
once during system design. The preamble based estimate of 

the channel matrix is therefore given by H = diag{li2}. 

Fig. QT| compares the performance loss of CP-OFDM and 
non-systematic coded UW-OFDM (G') in case the described 
preamble based channel estimate given by ( |67| i is used instead 
of perfect channel knowledge. As a highly interesting result 
we notice that both systems degrade by about the same scale: 
CP-OFDM experiences a loss of 0.8dB for r = 1 /2 and 0.6dB 
for r — 3/4, respectively (all results again measured at a BER 
of 1CP 6 ), while the performance of UW-OFDM degrades by 
0.7dB for both coding rates. 

VII. Conclusion 

In this work we expanded our recently introduced systematic 
coded UW-OFDM concept to non-systematic coded UW- 
OFDM. For that we introduced optimized code generator 
matrices that distribute the redundancy over all subcarriers 
instead of only over a dedicated set. We derived optimization 
criteria to find a class of code generator matrices that in 
case of AWGN conditions ensure minimum error variances 
on the subcarriers after the estimation process at the receiver. 
However, due to the advantageous combination of the channel 
coding and dispersive preprocessing abilities of the optimized 
code generator matrices, non-systematic coded UW-OFDM 



LU 



-* — CP-OFDM perfect 
-0 — CP-OFDM imperfect 
-* — UW-OFDM non-syst. coded perfect 
-e — UW-OFDM non-syst. coded imperfect 
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Fig. 1 1 . Impact of imperfect channel estimation on the BER performance 
of CP-OFDM and non-systematic coded UW-OFDM (C) in a frequency 
selective environment with additional outer coding (QPSK). 



particularly features its superior performance in frequency 
selective channels. We showed simulation results for selected 
code generator matrices in the AWGN case as well as in fre- 
quency selective environments. It turns out that non-systematic 
coded UW-OFDM impressively outperforms systematic coded 
UW-OFDM and classical CP-OFDM w.r.t. the spectral and the 
bit error ratio behavior. 
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Appendix A 
Derivation of the Global Minimum of Jblue 



In Appendix [A] we proof that the gradient of the 
cost function Jblue with respect to the vector s = 

1 T 

si S2 ••• SN d °f singular values of G is zero if and 
only if all singular values are identical. And we proof that 
every local minimum of the cost function Jblue is also a 
global minimum with JBLUE.min = ^d-^d/c- Using the SVD 
in ( tSTb we have 



tr{G H G} = tr{VS H U ff USV ff } = tr{VD 2 V ff } 

= tr{D 2 }, (68) 
(G^G)- 1 = (VS^U^USV^)- 1 = (VD 2 ^)- 1 



= V(D 2 )- 1 V H 



(69) 
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Inserting ((68) and d69l l into ( f30l > leads to the following 
expression for the cost function 

J BLVE = ^rtr{tf}tr{V(&)-iV H } 

cN d 



^Ltr{D 2 }tr{(D 2 )-*} 



cN, 




4 



(70) 



which can now be regarded as a function of the singular values 
of G. The gradient of Jblue with respect to s follows to 



ds 



cN d 



s N d 



-3 I „2 



S N 3 d ( s l 



'iV rf 



J 



(71) 



Setting the gradient to zero leads to the system of equations 



S N d 



2 i 2 

s x + ---s Nd . 



(72) 



It is easy to see that s\ = S2 = • • • sjv d := s solves 
the system of equations. Furthermore, by subtracting the i tb 
equation from the f h equation for all i,j G {1, ...,N d } 
with i 7^ j it immediately follows that si = s% = ■ ■ ■ SN d 
is in fact the only solution to this system of equations. 
Consequently, every possible candidate G for a local minimum 
satisfies G H G = s 2 I (cf. ( 1521 and its implications). Inserting 
G H G — s 2 I into the cost function d30T > leads to the same 
expression as in d48l that corresponds to the numerically found 
local minima. Hence, every G fulfilling G H G = s 2 I results 
in the same (and minimum) value JBLUE.min = v^Nd/c which 
therefore constitutes the global minimum of the cost function. 

Appendix B 

Derivation of the Global Minimum of Jlmmse 

In Appendix iBl we proof that 5Jlmmse/<?s = if and only 
if all singular values are identical. And we proof that every 
local minimum of the cost function Jlmmse is a ls° a global 
minimum with Jlmmse, min = G \^dl ( c + !)■ Inserting ( |68l and 
( |69l > into (|36T > leads to the following expression for the cost 
function 



Jlmmse — 0" d tr 



cN d 



tr{D 2 } 



VD 2 V ff 



(73) 



which can be regarded as a function of the singular values 
of G. Applying the matrix inversion lemma we immediately 
obtain 



-l 

2\-l \ -ktH 



Jlmmse - ^tr |l - V (v H V + ^^(D 2 )- 1 ) V 

(74) 

With this step we achieved, that every matrix to be inverted 
in d74T i has a diagonal structure. Having in mind that D 2 is a 
diagonal matrix with the squared singular values of G on its 
main diagonal, and after some rearrangements we obtain 



iV,, 



Jlmmse = a'd.Nd — a d cN d 



cN d s 2 + tv{B 2 }' 



(75) 



The partial derivation of the cost function Jlmmse with respect 
to the j th singular value follows to 

d Jlmmse 



dsj 



-2(j z d cN d s j x 



x E 



(cN d s 2 + tr{D 2 }) 2 {cN d s 2 + tr{D 2 }) 2 



(76) 



It is easy to see that d Jlmmse/ dsj — is fulfilled if Si = Sj 
for all i G {1, N d } with i ^ j. In fact si = ■ ■ ■ = SN d = s 
is the only solution to <9Jlmmse/<9s = 0. This can be proved 
by subtracting the equations resulting from d Jlmmse /dsi = 
and <9JLMMSE/<9sj = for all i ^ j, which is not difficult but 
a kind of exhausting. The remaining argumentation coincides 
with the one for the BLUE in Appendix |A] However, the 
expression for the global minimum Jlmmse. min = ^d-^W ( c +l) 
which is obtained by inserting G H G = s 2 I into (f36t differs 
from Jblue ,min- 
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